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O ■ ABSTRACT 
H ' 

\ We present the first INTEGRAL observations of tlie type 1.5 Seyfert galaxy NGC 4151. 

■ Combining several INTEGRAL observations performed during 2003, totaling ^ 400 ks of exposure 
time, allows us to study the spectrum in the 2 - 300 keV range. The measurements presented 
here reveal an overall spectrum from X-rays up to soft gamma-rays that can be described by an 

\ absorbed (iVn — 6.9 x lO^^cm"^) model based on a Compton continuum from a hot electron 

■ population [kT^ — 94keV) from an optically thick (r — 1.3) corona, reflected on cold material 
{R = 0.7), consistent with earlier claims. The time resolved analysis shows little variation of the 
spectral parameters over the duration of the INTEGRAL observations. The comparison with 
CGRO/OSSE data shows that the same spectral model can be applied over a time span of 15 
years, with flux variations of the order of a factor of 2 and changes in the underlying continuum 
reflected by the temperature of the electron population {kTe = 50 — 100 keV). When modeled 
with an exponentially cut-off power law plus Compton reflection this results in photon indices 
ranging from F = 1.5 to F = 1.9 and a cut-off energy in the range 100 — 500 keV. 

Subject headings: galaxies: active — galaxies: individual (NGC 4151) — gamma rays: observations — 
X-rays: galaxies — galaxies: Seyfert 
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1. Introduction 

The Seyfert galaxy NGC 4151 is one of the 
most extensively studied AGN (Active Galactic 
Nucleus). This is due to the fact that it is one 
of the brightest (V ~ 11 mag, bolomctric lumi- 
nosity ^ 10''"' erg s~^) and nearest (16.5 Mpc for 
Ho = 75kms^^) among its class. Observations 
throughout the electromagnetic spectrum revealed 
a complex and somewhat special case of a Seyfert 
1.5 galaxy (see e.g. Dermer & Gehrels (1995) and 
Ulrich (2000) for a general review). In the optical 
the host galaxy appears to be a large barred spi- 
ral with nearly face-on orientation, showing two 
thin spiral arms (Arp 1977). At soft X-rays NGC 
4151 Chandra data show a complex spectrum, in- 
cluding line emission from a highly ionized and 
variable absorber, in the form of H-like and He- 
like Mg, Si, and S lines, as well as lower ioniza- 
tion gas via the presence of inner-shell absorp- 
tion lines from lower-ionization species of these el- 
ements (Kraemer et al. 2005). Due to the heavy 
absorption, it is not possible to tightly constrain 
the intrinsic X-ray photon index below 10 keV. 
The harder continuum at 2 — lOkeV as studied 
by ASCA shows a flatter spectrum with F ~ 
1.5 (Weaver et al. 1994). Observations by XMM- 
Newton in the 2.5 — 12 keV energy range suggest 
that the spectrum can be modelled by an absorbed 
power-law continuum (fixed F = 1.65), a high en- 
ergy break fixed at 100 keV, a neutral Compton re- 
flection component with a reflection fraction of 2.0 
compared to the normalization of the power law, 
and a narrow iron Kq line. The fact that R > 1 
is explained by the relative variability timescales 
in the direct and reflected component. The ab- 
sorption of the power-law continuum was repre- 
sented by a product of two absorption components 
(Schurch et al. 2003). The Galactic hydrogen col- 
umn density in the line of sight is A^h — 2.1 x 
10^°cm~^. The spectrum detected by BeppoSAX 
seemed to be best described by a power law with 
an exponential cut-off at > 50keV and an iron 
fluorescence Ka line (Piro 1999; Piro et al. 2000). 
In addition to the fluorescence line, which re- 
veals the existence of cooler material surround- 
ing the AGN, a reflection hump in the 10 — 
30kcV has been reported in Ginga/OSSE data 
(Zdziarski et al. 1996), in combined ASCA and 
OSSE data (Zdziarski et al. 2002), and in Bep- 
poSAX data (Petrucci et al. 2001; Schurch & War- 



wick 2002). 

In this paper we present analysis of recent ob- 
servations of NGC 4151 performed by the INTE- 
GRAL (International Gamma-Ray Astrophysics 
Laboratory) satellite, and compare the results with 
previous studies. 

2. Simultaneous optical, X-ray and gamma- 
ray observations 

Observations in the X-ray to soft gamma-ray 
domain have been performed in late May 2003 by 
the INTEGRAL satellite (Winkler et al. 2003). 
This mission offers the unique possibility to per- 
form simultaneous observations over the 2 — 
8000 keV energy region. This is achieved by IN- 
TEGRAL'S Joint European X-ray Monitor (2-30 
keV) (JEM-X; Lund et al. 2003), the INTEGRAL 
Soft Gamma-Ray Imager (20-1000 keV) (ISGRI; 
Lebrun et al. 2003), and the Spectrometer aboard 
INTEGRAL (SPI; Vedrenne et al. 2003), which 
operates in the 20 - 8000 keV region. These three 
instruments use the coded aperture method for 
imaging. For a review on the particularities of 
this technique see Caroli et al. (1987). In ad- 
dition to these data an optical monitor (OMC, 
Mas-Hesse et al. 2003) provides V band photo- 
metric measurements. 

In order to compare the INTEGRAL results 
with another simultaneous mission, we use data 
from the Rossi X-ray Timing Explorer (RXTE). 

2.1. INTEGRAL data 

For the analysis of the high energy spectrum 
of NGC 4151 we used the data which were taken 
within 10 degrees radius around the source. Note 
that because of the nature of coded mask imag- 
ing, the whole sky image taken by the instrument 
has to be taken into account in the analysis, as all 
sources in the fleld of view contribute to the back- 
ground (Caroli et al. 1987). A journal of our ob- 
servations is shown in Table 1. The exposure are 
the ISGRI effective on-source times. This value 
is approximately the same for the spectrograph 
SPI, but the JEM-X and OMC monitors cover a 
much smaller sky area. Thus in the case of dither- 
ing observation, the source is not always in the 
field of view of the monitors. The analysis was 
performed using version 4.2 of the Offline Science 
Analysis (OSA) software distributed by the INTE- 
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GRAL Science Data Centre (ISDC; Courvoisier et 
al. 2003b). 

The analysis of the INTEGRAL/IBIS data is 
based on a cross-correlation proccdnrc between 
the recorded image on the detector plane and 
a decoding array derived from the mask pattern 
(Goldwurm et al. 2003). Imaging analysis of the 
ISGRI data finds NGC 4151 at R.A. = 12^10^32% 
decl. = 39°24'42" (J2000.0), 24" off the nominal 
optical position, with a significance of 224o-. The 
error radius of the position determination is 18". 
Note that this positioning accuracy is far below 
the instrumental angular resolution of 12' (see e.g. 
Gros et al. 2003). The error radius rules out the 
two other X-ray sources within 10', the LINER 
NGC 4156 and the BL Lac object MS 1207.9+3945 
(both at ~ 5' distance to NGC 4151 and both not 
detected by INTEGRAL'S instruments), to be the 
counter part of the hard X-ray emission. 

No cross-calibration correction had to be ap- 
plied between JEM-X and ISGRI, but the SPI 
data appear to be at a higher level by a factor of 
1.1. The SPI analysis was done using the specific 
analysis software (Diehl et al. 2003) including ver- 
sion 9.2 of the reconstruction software SPIROS 
(Skinner & Connell 2003) which is based on the 
"Iterative Removal of Sources" technique (IROS; 
Hammersley et al. 1992): a simple image of the 
field of view is made using a mapping technique 
which is optimised for finding a source assuming 
that the data can be explained by only that source 
plus background. The mapping gives the approx- 
imate location and intensity of the source, which 
are then improved by maximising a measure of 
the goodness of fit. The residuals of the fit are 
used as the input for a further image reconstruc- 
tion and source search (Skinner & Connell 2003). 
The SPI data further include some energy bins 
where SPIROS could not find a solution, and oth- 
ers which indicate problems in the background 
subtraction. This is not surprising as most of the 
data have been taken in staring mode, which is 
known to cause problems with background deter- 
mination in the SPI analysis (Dubath et al. 2005). 

The ISGRI spectra have been extracted using 
the standard method as described in Goldwurm et 
al. (2003): fiuxes and count spectra are extracted 
at the source positions for each pointing in prede- 
fined energy bins. Spectra of the same source col- 
lected in different pointings or "science windows" 



(SCWs) are then summed to obtain an average 
spectrum during the observation. Regarding the 
absolute flux calibration, it must be noted that 
the instrumental response of ISGRI is still not ac- 
curately determined. This problem is most severe 
for ofii'-axis observations. In addition the spectra of 
the Crab taken with JEM-X show systematic fea- 
tures in the 5 — 7 keV energy range, which might 
affect the measurement of the iron Ka line. 

With respect to the calibration unc;c;rtainty, the 
IBIS instrument team stated that the systematic 
error is of the order of 1.5% (2005, private commu- 
nication). Nevertheless, this value corresponds to 
on-axis observations within a short period of time 
with no disturbing influence, such as enhanced 
background activity. A combined fit of Crab spec- 
tra taken in revolution 43, 44, 120, 170, and 239, 
i.e. over a 1.5 year span of the mission, shows 
a larger uncertainty in the fiux. The scale of di- 
versity, assuming that Crab is a source with con- 
stant flux, gives some hint what is the scale of 
discrepancies in count rates observed in various 
conditions, i.e., with different dithering patterns 
or instrumental settings. Clearly the absolute cal- 
ibration cannot be more precise than the observed 
variations. Hence, it seems that for the 22-120 
keV band we may conclude that the joint la un- 
certainty of the spectral extraction method and 
the calibration files is about 5%. It may happen 
that for other sources and in some special circum- 
stances we may encounter a larger discrepancy, 
but for NGC 4151 observations after revolution 
65 and done mainly in staring mode a systematic 
error of 5% should be a valid approximation. 

The spectral shape calibration was also tested 
on INTEGRAL Crab observations. With the so- 
called canonical model for the Crab showing a 
single power law with photon index F = 2.1 
(Toor & Seward 1974). the values retrieved from 
INTEGRAL are P = 2.1 (.lEM-X; Kirsch et al. 
2004), and P = 2.2 for SPI and ISGRI. It should 
be kept in mind though that the Crab is signif- 
icantly brighter than NGC 4151, and systematic 
effects might depend on source brightness. 

NGC 4151 was too faint above ~ 150 keV to be 
detected by IBIS/PICsIT. Summing all available 

INTEGRAL data together shows that NGC 4151 
is detected up to 300 keV, in a total exposure of 
408 ks. 
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2.2. RXTE data 

The RXTE All Sky Monitor (ASM: Lcvinc ct 
al. 1996) scans about 80% of the sky every or- 
bit. This offers a unique way to monitor the emis- 
sion of bright X-ray sources like NGC 4151 in the 

I. 5 — 12keV energy range. We extracted fluxes 
from the RXTE/ ASM data base. The fltixes have 
been averaged over one day (see the triangles in 
Fig. 2). The data show a 10 ± SmCrab flux over 
the duration of the INTEGRAL observations, con- 
sistent with the JEM-X2 measurement. 

2.3. CGRO data 

The source was also monitored in the ~ 

20 - 100 keV band by the CGRO/BATSE ex- 
periment from 1991-2000 (Parsons et al. 1998; 
Harmon et al. 2004). Those measurements are 
made using an Earth-occultation technique, and 
were thus sampled one time per 90 min CGRO 
orbit. In practice, NGC 4151 is a faint source 
for that instrument; but daily summations lead to 
~ 3(7 determinations. We obtained the BATSE 
data spanning the ~ 3.5 yr overlap of the RXTE 
and CGRO missions, for the purpose of monitor- 
ing the long-term hard-to-soft spectral behavior. 

3. X-ray to Gamma-ray spectrum 
3.1. Time Averaged Spectral Analysis 

The overall spectrum extracted from the INTE- 
GRAL .JEM-X, ISGRl, and SPl data is shown in 
Fig. 1 together with the residuals of the model fit. 
A systematic error of 5% has been added to the 
ISGRI and SPl data. The model fitting results 
however do not depend strongly on the SPl data, 
because of the larger error bars compared to the 
ISGRI data. 

The spectral fitting was done using version 

II. 3.2 of XSPEC (Arnaud 1996). The data cover 
the energy range 2 — 300 keV. We applied a fit 
of an absorbed power law with an exponential 
high-energy cut-off, and added a Gaussian line 
in order to account for the iron Ka fluorescence 
line. The best fit values are F = 1.53+° °^, 
Ec = keV, and Nn = 5.1 +°;^ x 10^2 cm'^. 
The fluorescence line we obtain has an equiva- 
lent width of EW = 119 eV and a flux of /ko = 
4.5 +J g X IQ-^ phcm-2 (all errors are 90% con- 
fidence values). The statistics do not allow us 



to assess possible effects of rclativistic broaden- 
ing of the iron line (see e.g. Yaqoob et al. 1995). 
The source had an average flux of /20-200keV = 
0.012 phcm~2s~^ = 9.6 X 10^^"ergcm^2 s^^ dur- 
ing the observation. This model gives xt = 1-22 
for 196 degrees of freedom. The model of an ab- 
sorbed power law without a cut-off leads to a sig- 
nificantly poorer fit result {xl = 2.1). 

A physical model describing Comptonization 
of soft photons by a hot plasma, the so called 
compTT model, has been developed by Titarchuk 
(1994). The model includes the plasma tempera- 
ture Te of the hot corona, the optical depth t of 
this plasma, and the temperature Tq of the soft 
photon spectrum. Because the spectrum starts at 
2 keV, To is not well constrained by the data and 
has been fixed to 10 eV. The model fit resulted in 
Nn = 6.4+°;^ x 10^2 cm'^, kT^ = 39+" keV, and 

r = 1.37+o;33 for plane geometry. The compTT 

model gives xt = 1-29 for 195 degrees of freedom. 

Because those models did not give conclusive fit 
results, we added a reflection component from cold 
material to the cut-off power-law (the so-called 
PEXRAV model; Magdziarz & Zdziarski 1995). 
Prom radio measurements the inclination angle i 
of the accretion disk is estimated to be i ~ 65° 
(Pedlar et al. 1993), where i = would be the 
case of a face-on disk. The free parameters in 
the PEXRAV model are best fitted with a pho- 
ton index of F = 1.85 +g Qg, a hydrogen column 
density of TVh = 6.9 to^^ x 10^^ cm'^, a cut-off en- 
ergy of Ec = 447+^gQ keV and a relative reflection 
of i? = 1.0 +"'3. The large uncertainty and high 
value of the cut-off energy, which is well outside 
the energy range covered by our detection, indi- 
cates that this parameter is not well constrained. 
The PEXRAV model leads to xl = 0.99 for 195 
degrees of freedom. 

Although this model is a good representation of 
th{^ data, the choice of a power law for the underly- 
ing high-energy continuum is purely phenomeno- 
logical. A more physically approach involves 
the compPS model (Poutanen & Svensson 1996). 
This is a self-consistent model which computes 
a Comptonization spectra using exact numeri- 
cal solutions of the radiative transfer equation. 
The resulting spectrum is reflected from the 
cool medium according to the computational 
method of Magdziarz & Zdziarski (1995). The 
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model fit resulted in an electron temperature 
of kTe = 94^j*p keV, an optical depth of the 
corona r = 1.31 ^q q^, and a relative reflection 
ofR = 0.72 with xl = 0.98 for 196 degrees 
of freedom. 

We also applied models used in previous stud- 
ies. For example, using the model described in 
Schurch ct al. (2003) with the continuum fixed to 
an absorbed power law with F = 1.65, the fit re- 
sult worsens to xl = 1-06, and when also fixing 
the high energy break to 100 keV the fit results in 
xi = 1.39. It should be pointed out that though 
Schurch et al. used ASCA, BeppoSAX, and XMM- 
Newton data to develop their spectral 'template' 
model, and these data do not constrain a high- 
energy cut-off above 100 keV particularly well. 

Using different strength of the reflection param- 
eter, as described in Petrucci et al. (2001) also led 
to poorer fit results. For R = [0.0,0.01,0.2,0.45] 
the fit resulted in xl = [1.22,1.21,1.11,1.03]. A 
reflection component with R> 0.5 is therefore re- 
quired to achieve an acceptable fit result to the 
INTEGRAL data. 

Although previous studies showed that the spa- 
tial and dynamical structure of the absorbing me- 
dia in NGC 4151 is more complex (Kraemer et al. 
2005; Ogle et al. 2000; Schurch & Warwick 2002; 
Schurch et al. 2003) than the single cold absorber 
applied in our model, and that the soft X-rays are 
dominated by radiative recombination continua 
and X-ray emission lines (Kraemer et al. 2005; 
Schurch et al. 2004), a more complicated model 
incorporating such would not improve the fit re- 
sults to the INTEGRAL data. This is due to the 
fact that the JEM-X data in the relevant range up 
to ^ 3 keV do not have sufficiently high signal-to- 
noise to derive constraining information from this 
part of the spectrum. 

3.2. Spectral Evolution 

In order to study the evolution of the spec- 
trum in time, we applied ten time bins over the 
total observation period, combining about 33 ks 
of effective on-time in each bin. We used only 
those data for which we have simultaneous JEM- 
X and ISGRI coverage. We then applied the cut- 
off power-law model to the individual time bins 
for ISGRI and an absorbed power law model for 
JEM-X and extracted fluxes in the 20 - 100 keV 



and 1.5 — 12kcV range, respectively. The results 
are shown in Fig. 2 together with the RXTE/ ASM 
lightcurve in the same energy band as JEM-X. 
No significant variation of the flux can be seen 
here, and also the fit parameters are consistent 
over the observed period within the error bars. 
In order to look for variability on shorter time 
scales we then binned the INTEGRAL JEM-X2 
(5 - 20 keV) and ISGRI (20 - 100 keV) data into 
2000 s time intervals (Fig. 3). Both lightcurves 
show an apparent flare at f ~ 1243.0 IJD^. A 
more careful analysis of the data enclosed within 
that time bin, 1242.5 IJD < t < 1243.5 IJD, 
does show a slightly increased flux (/20-200keV = 
0.0132 phcm-^s"^ = 1.03 x IQ-^ergcm-^s"^) 
but the spectral parameters are consistent with 
the time-averaged values. 

We performed a cross-corrleation analysis of the 
JEM-X and ISGRI light curves (Fig. 3), apply- 
ing the "Discrete Cross-Correlation" algorithm of 
Edelson & Krolik (1988). A lag-spectrum span- 
ning ±0.5 days, with a granularity of 0.05 days 
(or roughly 2.5 times our sampling time) was con- 
structed. A plot of the resulting correlation statis- 
tic versus lag was symmetric about 0, with a peak 
value of about 40 %. This suggests that the hard 
and soft X-ray variations are well correlated, but 
that our sampling timescales are too coarse to 
study any reprocessing timescale on or light travel- 
time effects. 

We also evaluated the hard-X-ray spectral vari- 
ability following the procedure of Vaughan et 
al. (2003) and Markowitz, Edelson, & Vaughan 
(2003). We computed the fractional rms variabil- 
ity amplitude based on the per-SCW coimt rates 
and uncertainties on each of 4 energy bands: 13- 
20, 20-36, 36-58, and 58-87 keV. Our results are 
limited by the data: there were a total of 244 
SCWs spanning about 7 days. The signal-to-noise 
ratios ranged from ~ 17 in the best case (20-36 
keV) to ^ 3 — 5 in the poorest case (13-20 kcV). 
For purpose of comparison, we have reproduced 
the lower energy results of Markowitz, Edelson, & 
Vaughan (2003), and presented with our results 
in Fig. 4. 

For the same time intervals as used in the 
high-en(-r,2,y analysis, optical photometry values 

^LJD is related to the modified Julian Date by IJD = 
MJD - 51544.0 
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from the OMC have been extracted. Due to 
the high signal-to-noise these data have been 
binned into 1000 s time intervals (Fig. 5). We 
also considered possible correlation between vari- 
ations in our optical measurements and the ISGRI 
hard X-ray rates. We derived the approximate 
mono-chromtatic flux (f °) from the OMC 

^■'5500 A' 

V-band magnitudes, and normalized these val- 
ues to 10~^^ ergscm"^ s~^ A~^. We then apphed 
the same cross-correlation algorithm to the X-ray 
count rates and the optical fltix. The same range 
and granularity as we used for the JEM-X - ISGRI 
case was applied. The results obtained indicate 
a lack of significant correlation (peak amplitude 
^ 25%), and a flat, symmetric distribution about 
zero-lag. 

3.3. Search for 7-ray line emission 

The SPI spectrometer offers the opportunity 
to look for narrow line features. In case the 
high energy spectrum would be dominated by 
non-thermal processes, one could expect a sig- 
nificant amount of e^ pair annihilation processes 
around 511 keV. At this energy, SPI has a reso- 
lution of FWHM = 1.95 keV (Attie et al. 2003). 
CGRO/OSSE was able to determine the 3cr up- 
per limit of a broad (450 — 600 keV) emission line 
feature with / < 6 x 10~^ phcm~^ s~^. We there- 
fore looked for a comparably narrow feature by 
extracting an image in the 499 — 519 keV en- 
ergy band, centred on the 511 keV rest-frame 
energy at NGC 415rs rcdshift. No emission 
at the restframe energy of the annihilation line 
is detectable. The 3a upper limit extracted 
from the SPI data by SPIROS is .Aog-siofeeV = 
lO"** phcm~^ s~^. The 3a line sensitivity derived 
from ground calibration gives a similar value of 
6.5 X 10-'"^phcm-2s-i (Roques et al. 2003). The 
measurement presented here is therefore not im- 
proving the CGRO/OSSE upper limit, because of 
the much shorter INTEGRAL/SPl exposure time. 

4. Discussion 

The overall INTEGRAL spectrum can be mod- 
elled by Comptonization of photons in a hot 
plasma, and suggest that such an environment 
characterizes the central regions of NGC 4151. 

The best comparable coverage in the energy 
range studied by INTEGRAL ISGRI and SPI has 



been provided by CGRO/OSSE (Johnson et al. 1997). 
The measurements from INTEGRAL and OSSE 
are compared in Tab. 2. The OSSE results have 
been derived by fitting simultaneously the individ- 
ual OSSE spectra, allowing for different normali- 
sation. Our results are similar to what has been 
reported in Johnson et al. (1997). We compared 
the OSSE data directly with the INTEGRAL ob- 
servation through simultaneous fitting. Figure 6 
shows a simultaneous fit of both data sets, the 
average INTEGRAL and average CGRO/OSSE 
spectrum. The data can be represented by a 
(non-physical) cutoff power law with individual 
constant factors in order to account for fiux vari- 
ability (xt = 1-10). Note that the energy range 
above ^ 100 keV is dominated by the OSSE data, 
as their error bars are significantly smaller than 
the INTEGRAL ones, because the INTEGRAL 
exposure time is much smaller. The fit shows that 
the 20—200 keV flux varied by a factor of 2 over the 
individual INTEGRAL and GGRO observations. 
It is remarkable that the same simple model with 
F = 1.5 and Ec = 95keV can represent the data 
from the different CGRO/OSSE and INTEGRAL 
observations. 

When applying the more meaningful PEXRAV 
or compPS model to the CGRO data it turns out 
that a reflection component is not constrained by 
the OSSE data. This is not surprising, as there 
are no OSSE data below 50 keV where a reflection 
component would have a signiflcant influence. The 
PEXRAV flt to the combined data shows that a 
reflection component improves the fit according to 
an F-test. The reduced is 1-07 and the relative 
reflection is only R = OA'^^''^ in this case. The 
compPS model gives a better flt result (x^ = 1-01 
for 476 d.o.f.) with an electron temperature of 
kTe = 55 '^'^ keV and a relative reflection of i? = 

-0.1- 

These results are mainly consistent with previ- 
o\is BeppoSAX observations. Three observations 
by BeppoSAX, two in 1996 and one in 1999, sug- 
gest that the strength of the reflection component 
is variable, with -Riggeif = 0.01,i?i999 — 0.2, and 
-R1996L = 0.45 (Petrucci et al. 2001). R is nor- 
malized so that ii = 1 in the case of an isotropic 
source above an infinite reflecting plane. R is of- 
ten considered as an estimate of ^/2n where fl 
is the solid angle subtented by the reflector as 
seen from the isotropic X-ray source. It has to be 
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pointed out that a relative reflection of i? > 1, 
as seen in the work of Schurch et al. (2003), 
represents the case when more primary X-ray ra- 
diation is emitted toward the reflector than to- 
ward the observer. This can be explained by vari- 
able emission of the central engine and a time 
delay between the underlying continuum and the 
reflected component, caused by a large distance 
of the reflecting material to the primary source 
(Malzac & Petrucci 2002). Another model which 
can lead to large R values is based on relativis- 
tic eff'ects caused by a dynamic corona moving 
towards the reflecting disk (Malzac ct al. 2001). 
Other explanations include a special geometry 
with a high intrinsic covering fraction of the cold 
disk material (Malzac 2001) and general rclativis- 
tic light bending effects caused by a maximally ro- 
tating Kerr black hole (Miniutti & Fabian 2004). 

Another study of the high energy spectrum of 
NGC 4151 combining CGRO/OSSE and ASCA 
data has been presented by Zdziarski et al. (2002) . 
This work came to the conclusion that the data 
can be represented by a simple cut-off power 
law with a narrow iron Ka line, adding ther- 
mal bremsstrahlung at low energies, which are 
not covered by our data. This flt resulted in a 
= 1-6 and the authors added a broad Gaus- 
sian line in order to account for residuals in the 
5 — 6.5 keV energy band. Zdziarski et al. argue 
that a model with complex absorption and Comp- 
ton reflection component is more likely to give a 
physical meaningful explanation for the observed 
spectrum. They also discovered that the absorber 
in NGC 4151 shows evidence for the presence of 
several absorbing components in the line of sight, 
which cannot be conflrmed by the INTEGRAL ob- 
servations because of the lack of data below 2 keV. 

Matt ct al. (2003) argue that absorption in 
Seyfert 2 galaxies can be caused by several circum- 
nuclear regions. They favour a model with a pos- 
sible temporary switching-off of the nuclear radi- 
ation to explain observed variability. Their model 
predicts an iron Ka equivalent width of ~ 100 cV 
for a Compton thick case with A^h — 10^"' cm^^, 
consistent with the measurements for NGC 4151. 
This could mean that the iron line is produced in 
material in the line of sight, at large distances from 
the central engine. 

Schurch & Warwick (2002) showed that the 
same 'template' spectrum matches observations 



of NGC 4151 done by ASCA, BeppoSAX and 
even XMM- Newton (Schurch et al. 2003). Their 
model spectrum includes a strong reflection com- 
ponent (i? ~ 2 from PEXRAV) as well as an ab- 
sorption component with variable iron abundance. 
The INTEGRAL data are not suitable for detailed 
study of the absorption components. The obser- 
vations show that the underlying continuum is 
significantly steeper, with F = 1.85 ^Q ^g instead 
of r = 1.65+°;°3 (Schurch & Warwick 2002), and 
with a higher cut-off energy, although the flt to 
the INTEGRAL data using F = 1.65 is still ac- 
ceptable {xi = 1-06 instead of 0.99). This shows 
that it is possible to represent the spectrum of 
NGC 4151 over a long time scale with the same 
model with only a few parameters showing signifl- 
cant changes. The PEXRAV model appears to be 
a good representation of the high energy spectrum 
of NGC 4151, with the continuum's photon index 
varying in the range r = 1.5. ..1.9, and the cut-off 
energy is Ec > 90keV with no clear upper limit 
detected at present. 

The variability study we performed on the IN- 
TEGRAL data showed that the hard and soft X- 
ray variations are well correlated during the ob- 
servation. The fractional rms variability ampli- 
tude based on the per-SCW count rates and un- 
certainties in four energy bands in comparison to 
an earlier study by Markowitz et al. (2003) sug- 
gests that there is a continuation of the weak anti- 
correlation, or even a flattening, of the variabil- 
ity with energy (Fig. 4) .While one must bear in 
mind that these data are from separate epochs, 
they are consistent with a scenario where the puta- 
tive Compton-reflected component is slowly vary- 
ing relative to the disk emission. This in turn is 
suggestive of a reflector which may be larger in 
scale and/or farther from the low-energy photon 
sources than the inner accretion disk; perhaps as- 
sociated with an outer disk torus, or with outflow- 
ing plasma. 

Most Seyfert galaxies show a softening of the 
X-ray continuum as sources brighten (Markowitz, 
Edelson, & Vaughan 2003; Nandra et al. 1997). 
We did not see signiflcant variations of the spectral 
parameters over the duration of our observations. 
In particular, when modelling the spectrum with 
a cut-off power-law, we did not detect any cor- 
relation of flux with spectral slope. This can be 
understood when studying the correlation of spec- 
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tral hardness ratio with flux as derived from simul- 
taneous CGRO/BATSE and RXTE/ ASM data as 
shown in Fig. 7. Here we have plotted the hardness 
ratio (20 - 200 keV) / (2 - 10 keV) versus the RXTE 
(2 — 12keV) flux based on daily averages from 
1996-2000. An overall correlation is discernible, 
but the scatter is large and thus might not be de- 
tectable over a time span of 400 ks (as in our case), 
and can even be missed when comparing the long- 
term observations by CGRO/OSSE. 

5. Conclusion 

The results from the INTEGRAL JEM-X2, IS- 
GRI, and SPI data show that the high-energy 
spectrum of NGC 4151 can be approximated by 
a rather simple model including Compton reflec- 
tion. A Compton reflection model (PEXRAV) 
results in an underlying power-law with F = 
1.85 ^Q Qg, a hydrogen column density of A^h = 
6.9+0 4 X 10^2 cm-2, a cut-off energy of Ec = 
450+200 keV, ^'^^ ^ relative reflection of i? = 
1-0 +0 3 • However, the high-energy cut-off is still 
not well constrained by the data. The more phys- 
ical compPS model reveals a hot electron popula- 
tion (kTe = 94 keV) , an optically thick corona 
(t = 1.3+"'^), and a significant reflection com- 
ponent {R = 0. ?+[] :[). Even though the data do 
not constrain the reflection model strongly, they 
show the necessity of a reflection component. The 
study shows that not only the reflection compo- 
nent varies on long time scales, but that also the 
underlying continuum changes signiflcantly over 
the years. 

The comparison to previous CGRO/ OSSE mea- 
surements shows that this model is a valid repre- 
sentation of the high energy spectrum over decadal 
timescales, with the normalisation varying by a 
factor of 2 and the temperature of the underlying 
continuum varying in the range kTe = 50—100 keV 
for the various observations. A significant varia- 
tion of the spectral parameters during the INTE- 
GRAL observation cannot be detected. The ef- 
fect of spectral softening with higher flux values 
can only be seen in observations covering longer 
time scales. This is shown in the comparison of 
CGi? 0/BATSE (20-200 keV) to RXTE/ ASM (2- 
10 keV) fluxes. 

The INTEGRAL data cannot confirm the ex- 



istence of a high iron overabundance in the ab- 
sorber, simply because of the lack of data below 
2keV. 

The thermal Comptonization is seen in the cut- 
off of the spectrum, ruling out synchrotron emis- 
sion being dominant in this source. This is also 
supported by the non-detection of emission in the 
annihilation line in NGC 4151. 

The study of NGC 4151 indeed shows that the 
INTEGRAL mission can give valuable insights on 
the characteristics of the high-energy cut-off in 
bright AGN, as has been already shown for the 
Seyfert 2 NGC 4388 (Beckmann et al. 2004a, 
2004b). This shows that AGN can be weU de- 
scribed by an optical thick disk with a hot corona. 
We will present a general overview on the proper- 
ties of AGN observed by INTEGRAL in a forth- 
coming paper (Beckmann et al. 2005). 
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NGC 4151 INTEGRAL spectrur 

compPS model 
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Fig. 1.— Summed spectrum of the INTEGRAL 
JEM-X2, ISGRI, and SPI data {upper panel). The 
data are well represented by the Compton reflec- 
tion model compPS with electron temperature of 
kTe — 94keV and relative reflection of i? = 0.7 . 
In addition a narrow gaussian iron Ka line with 
equivalent width EW = 119 eV appears in the 
spectrum. The lower panel shows the residuals 
of the fit. 



This 2-column preprint was prepared with the AAS lAT^rjX 
macros v5.2. 
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Fig. 2.— INTEGRAL and RXTE/ ASM light 
curves in Crab units. The squares and octagons 
represent the 20 - 100 keV ISGRI and 1.5-12 keV 
JEM-X2 hghtcurve, respectively. Fluxes were ex- 
tracted from model fits to the data. The triangles 
are RXTE/ ASM 1.5 - 12keV daily averages. 
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Fig. 4. — Fractional rms variability amplitude 
based on the INTEGRAL ISGRI count rates {cir- 
cles) compared to the results of Markowitz, Edel- 
son, & Vaughan (2003; triangles). 
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Fig. 3.— INTEGRAL ISGRI {upper curve, 2000 s 
binning) and JEM-X2 (1000 s binning) lightcurve 
in the 20 - 100 keV and 5-20 keV energy band, 
respectively. 



11.1 - 




1240 1242 1244 

time [IJD] 



Fig. 5.— INTEGRAL/OMC light curve. Error 
bars are la values. 
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Table 1: INTEGRAL observations 



obs. date 


INTEGRAL 


exposure 


number of 


observation 


start 


revolution 


time [ks] 


pointings 


mode 


22/05/03 


73 


7.3 


3 


5x5 


23/05/03 


74 


170.8 


10 


5x5, staring 


25/05/03 


75 


174.5 


5 


staring 


28/05/03 


76 


55.1 


4 


staring 



Table 2: Spectral fits to average INTEGRAL (2 - 300 keV) and OSSE (50 - 600 keV) spectra 



Model 


Mission 


Photon 

Index 


Ec or kT 
[keV] 


T or 


Xtidof) 


Power law 


INTEGRAL 


1.79 ±0.02 






2.10 (197) 


Power law 


CGRO/OSSE 


9 +0.03 
^•00 _o.03 






1.73 (278) 


PL-exp 


INTEGRAL 


1 rq +0.04 
_o.04 


124 til 




1.22 (196) 


PL-exp 


CGRO/OSSE 


1 47 +0-15 

-0.16 


90 t\l 




0.98 (277) 


PL-exp 


INTEGRAL+CGRO 


1 49 +°°^ 

'-■^^ -0.04 


94 tl 




1.10 (476) 


compTT" 


INTEGRAL 




39 +!r 


1 07 +0.26 
-0.33 


1.29 (195) 


compTT" 


CGRO/OSSE 




44 +6 


89 +°-^^ 

-0.17 


0.99 (276) 


PEXRAV 


INTEGRAL 


-0.09 


447 +885 

-190 


98 ^"-^"^ 

-0.30 


0.99 (195) 


PEXRAV 


CGRO/OSSE 


1 rc9 +0.25 
-0.28 


94 +^9" 


00 ^^■''^ 


0.99 (276) 


PEXRAV 


INTEGRAL+CGRO 


1 61 +°°^ 
'-■^^ -0.06 


122 til 


40 +°-^^ 

I'-'*" -0.08 


1.07 (475) 


compPS'' 


INTEGRAL 




94 +^ 

-10 


R = 0-72tllt 


0.98 (196) 


compPS'' 


CGRO/OSSE 




52 +1 


i? = 0.48+°-^^ 


0.98 (277) 


compPS*" 


INTEGRAL+CGRO 




55 +^ 


R = OM+_ilt 


1.01 (476) 



" Following Titarchuk 1994 

^ Following Poutanen & Svensson 1995 



12 





10 100 
chonnel energy (keV) 

Fig. 6.— INTEGRAL and CGi?0/OSSE simul- 
taneous data fit. The OSSE data represent the 
average of the NGC 4151 spectra. All data are 
well represented by a Compton reflection model 
(compPS; Poutanen & Svensson 1996), with the in- 
dividual OSSE data scaling with 0.6 . . . 1.2 versus 
the INTEGRAL/ 3'EM-X2 data. 




Fig. 7. — Spectral hardness derived from the 
BATSE and RXTE data, versus the RXTE/ ASM 
flux. A softer spectrum is correlated with a higher 
flux. The error bars on the lower left show the 
typical uncertainties of the measurements. 
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